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Abstract. The gravitational lensing of gravitational waves should be treated in the wave optics instead of the geometrical optics 
when the wave length A of the gravitational waves is larger than the Schwarzschild radius of the lens mass M. The wave optics is 
based on the diffraction integral which represents the amplification of the wave amplitude by lensing. We study the asymptotic 
expansion of the diffraction integral in the powers of the wave length A. The first term, arising from the short wavelength limit 
\ /I — * 0, corresponds to the geometrical optics limit. The second term, being of the order of A/M, is the leading correction term 

^ . arising from the diffraction effect. By analysing this correction term, we find that (1) the lensing magnification /i is modified to 

Ijr^ ' (1 + S), where S is of the order of (A/M)-, and (2) if the lens has cuspy (or singular) density profile at the center p(r) oc r"" 

, (0 < Q- < 2), the diffracted image is formed at the lens center with the magnification // ~ (A/My^". 
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1. Introduction 



(N 
O 

o 

^p^, The gravitational lensing of light is usually treated in the geometrical optics approximation, which is valid in almost all obser- 
O ' vational situations since the wave length of light is much smaller than typical scales of astrophysical lens objects. However for 
^ the gravitational lensing of gravitational waves, the wavelength is long so that the geometrical optics approximation is not valid 
^ in some cases. As shown by several authors (Ohanian 1974, Bliokh & Minakov 1975, Bontz & Haugan 1981, Thorne 1983, 
• • , Deguchi & Watson 1986), if the wavelength A is larger than the Schwarzschild radius of the lens mass M, the diffraction effect 
becomes important. Thus, the diffraction effect is important for the lens mass smaller than 10** M© (/I/IAU), where 1 AU is the 
, wavelength for the planed laser interferometer space antenna (LISA: Bender et al. .2000. 1 
H I From the above discussion, for M the diffraction effect is important and the magnification is small (the wavelength is so 
■ - - long that the wave does not feel the existence of the lens), and for /I <K M the geometrical optics approximation is valid. In this 
paper, we consider the case for A ^ M, i.e. the quasi-geometrical optics approximation which is the geometrical optics including 
corrections arising from the effects of the finite wavelength. We can obtain these coiTection terms by an asymptotic expansion of 
the diffraction integral in powers of the wavelength A. ' The diffraction integral represents the amplification of the wave amplitude 
by lensing in the wave optics. It is important to derive the correction terms for the following two reasons: (1) calculations in the 
wave optics are based on the diffraction integral, but it is time consuming to numerically calculate this integral especially for high 
frequency (see e.g. Ulmer & Goodman 1995). Hence, it is a great saving of computing time to use the analytical expressions. (2) 
We can understand clearly the difference between the wave optics and the geometrical optics (i.e. the diffraction effect). 

This paper is organized as follows: In §2 we briefly discuss the wave optics in gravitational lensing of gravitational waves. In 
§3 we show that in the short wavelength limit /I — > 0, the wave optics is reduced to the geometrical optics limit. In §4 we expand 
the diffraction integral in powers of the wavelength A, and derive the leading correction terms arising from the effect of the finite 
wavelength. In §5 we apply the quasi-geometrical optics approximation to the case of the simple lens models (the point mass 
lens, SIS lens, isothermal sphere with a finite core lens, and the NFW lens). Section 6 is devoted to summary and discussions. 
We use units of c = G = 1 . 



Send offprint requests to: R. Takahashi 

' The asymptotic expansion of the diffraction integral has been studied in optics. See the following text books for detailed discussion: Kline 
& Kay (1965), Ch.XII; Mandel & Wolf (1995), Ch.3.3; Born & Wolf (1997), App.III, and references therein. 
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2. Wave Optics in Gravitational Lensing 

We consider gravitational waves propagating near a lens object under the thin lens approximation in which the gravitational wave 
is scattered only on the thin lens plane. The gravitational wave amplitude is magnified by the amplification factor F which is 
given by the diffraction integral as (Schneider et al. I1992> . 

F(u, 71) = ^ f d^^ exp [icotdi^, I/)] , (1) 

where oj is the frequency of the wave, ^ is the impact parameter in the lens plane, ij is the source position in the source plane. 
Dl, Ds and Dls are the distances to the lens, the source and from the source to the lens, respectively. The time delay from the 
source position i/ through ^ is given by, 

td(^^n)-^^lTr-4-] -kO + kriri). (2) 



2Dls \Dl Ds, 

The deflection potential 0'(^) is determined by, V^i^ = SttS, where SC^) is the surface mass density of the lens. 0,„ is the additional 
phase in F, and we choose 0,„ so that the minimum value of the time delay is zero. F is normalized such that \F\ = 1 in no lens 
limit ((A = 0). 

It is useful to rewrite the amplification factor F in terms of dimensionless quantities. We introduce as the normaUzation 
constant of the length in the lens plane. The impact parameter ^, the source position r/, the frequency and the time delay are 
rewritten in dimensionless form, 

€ Dl Ds 2 rr DiDis 2. z-on 

so goi^s DlsUi Us 

We use the Einstein radius (~ (MD)'^^) as the arbitrary scale length ^0 for convenience. Then, the dimensionless frequency is 
w ~ Ml A from Eq.0. The dimensionless time delay is rewritten from Eq.Q and Q as, 

r(x,y)= i|x-y|2-,A(x) + 0„,(y), (4) 

where iA(x) and 0m(y) correspond to and 0„,(77) in Eq.©: (lA, 4'm) - DiDis /{Ds^q) (iA, 4>m)- Using the above dimensionless 
quantities, the amplification factor is rewritten as, 

Fiw, y) = ^ J" exp [iwT(x, y)] . (5) 
For the axially symmetric lens models, the deflection potential iA(x) depends only on x = |x| and F is rewritten as. 



f 

Jo 



F(w,y) - -iwe'"-'^^^ I dx x Joiwxy) exp 



iwl ]-x^ - i//{x) + 4>„,(y) 



(6) 



where Jq is the Bessel function of zeroth order It takes a long time to calculate F numerically in Eqs.(|5} and (|6j, because the 
integrand is rapidly oscillating function especially for large w. In the Appendix A, we present a method for numerical computation 
to shorten the computing time. 



3. Geometrical Optics Approximation 

In the geometrical optics limit (w » 1), the stationary points of the r(x, y) contribute to the integral of Eq.Q so that the image 
positions Xj are determined by the lens equation, V^T(x, y) = or y = x - V^^/^Cx). This is just the Fermat's principle. We expand 
T(x, y) around the j-th image position x^ as, 

T(x, y) = r,- + i 2 dad.Tixj, y)x„x, + O(x^) (7) 

^ a,b 

where x = x - x^, Tj - T{Xj, y), and the indices a,b,... run from 1 to 2. Inserting Eq.0 to (|5|i, we obtain the amplification factor 
in the geometrical optics limit as (Schneider et a/. 1 19921 Nakamura & Deguchi ll999l l. 

F,Uw,y)^J]\^'Jf" ^'"'^'''""'^ (8) 
j 

where the magnification of the j-th image is jij - Ij det(5y/(9xy) and rij -0,1 /2, 1 when Xj is a minimum, saddle, maximum 
point of r(x, y). 



R. Takahashi: Quasi-geometrical Optics Approximation in Gravitational Lensing 



3 



4. Quasi-geometrical Optics Approximation 

4.1. Effect on the Magnification of tlie Image 

We expand the amplification factor F(w, y) in powers of 1/w (<K 1) and discuss the behavior of the order of 1/w term. Here, we 
only consider axially symmetric lens models because the basic formulae are relatively simple, while the case of the non-axially 
symmetric lens models is discussed in the Appendix B. We expand T(x, y) in Eq.0 up to the fourth order of x as. 



r(x,y) = Tj + ajx\ +l}jx\ + ^ ^ dadbdcT(Xj,y)xaXbXc + ^ ^ dadhdcddT(xj,y)xaXbXcXd + O(x^), 

a,h,c a,b,c,d 

where aj and j3j are defined by, 
a, = i(l-^;), = i 

where ip^p - d"il/(\Xj\)/dx". Inserting Eq.(|5) into (|5) and expanding Eq.(|5) in powers of 1/w, we obtain, 
^^^'^^ = iS^'""' r l + -i-Y^dMT(xj,y)x:xlx: + U~^ 

J L ' \a,b,c 

+ 4 Z dadhdATixj,y)x'^x',x'A+0(w-^'^) 



(9) 



(10) 



a.b.c.d 



(11) 



where x' - y/wx. The first term of the above equation il It is the amplification factor in the geometrical optics limit Fgeo in 
Eq.ljSj. The integral in the second term vanishes because the integrand is an odd function of x'^. The third term is the leading 
correction term, being proportional to l/w, arising from the diffraction effect. Thus the deviation from the geometrical optics is 
of the order of l/w ~ A/M. Inserting T(x, y) in Eq.@ into we obtain F as. 



(12) 



where 



A; = — 

■' 16 



1 



(3) 

2 1 -A) , aj-/3j I 



afij \xj 

and Aj is a real number. We denote dF,„ as the second term of Eq. (ll2t . 

dF„,(^v,y)^-y^j\^ij\''^e"'^'-i'''y. 



(13) 



(14) 



Since dF,n is the correction term arising near the image positions, this term represents the corrections to the properties of 
these images such as its magnifications, and the time delays. We rewrite F in above equation Mil as, F{w,y) - 2y| yU;[l + 
(Aj/w)^]|'^^e"^'-^'^'*'^j"'™j + (9(w"^), where 6ifj - arctan(Aj/w). Thus in the quasi-geometrical optics approximation, the magni- 
fication Hj is modified to //^[l + (Aj/w)^], where (A^/w)^ is of the order of {A/My^. That is, the magnification is slightly larger 
than that in the geometrical optics limit. The phase is also changed by 6(pj, which is of the order of A/M. 



4.2. Contributions from the Non-stationary Points 



In the previous section, we showed that the contributions to the diffraction integral F arise from the stationary points (or image 
positions). In this section, we discuss the contributions from the non-stationary points. We denote x„s as the non-stationary point, 
at which the condition |V vTI ^ satisfies. If T{x, y) can be expanded at x„s, we obtain the series of T similar to Eq.0 as. 



r(x, y) = T„s + T[xi + T2X2 + a„sX^ + Pnsx\ + 0{x^), 



(15) 



where x = x - x„s, T„s - T(x„s, y) and 2 ~ t^i,27"(x,,j, y). Note that either T[ or does not vanish because of |Vr| at x„ 
a„s, and /3„s are defined same as Ea.( llO> . but at the non-stationary point x„s. 

Inserting Ea.( ll5> to l|5|l and expanding the integrand in powers of l/w, we obtain. 



F(w,y) — 2n- 



iw 



6(T[)6m) - - 
w 



ans 



02 \ 



6{T[)6(r2)+0{w-^) 



(16) 
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Thus, if both TJ + and + 0, the above equation vanishes. Even if either = or - 0, it is easy to show that F vanishes. 
Thus, the contributions to the amplification factor F at the non-stationary points are negligible. 

In the above discussion, we assume that r(x, y) has the derivatives at the non-stationary point x„j in Ea.(ll5>. But, if the 
derivatives of T are not defined at x„j, the result in Ea. (ll6> should be reconsidered. If the lens has the cuspy (or singular) density 
profile at the center, the derivatives of T are not defined at the lens center. We will discuss this case in the next section. 

4.3. Central Cusp of the Lens 

We consider the correction terms in the amplification factor F arising at the central cusp of the lens. For the inner density profile 
of the lens p cc r " (0 < a < 2),^ the surface density and the deflection potential at small radius are given by, 

2(f) oc ^-"^^ forff 1, 

oc In^ for a - I, (17) 
i/r(x) oc x-"*^ fora 5t 1, 

oc x^lnx for a - 1. (18) 

We note that the Taylor series of i/'(x) around x = like Ea.( ll5> cannot be obtained from the above equation ( I18l l. For example, 
in the case of a = 2, the deflection potential is oc x, but the derivative of |x| is discontinuous at x = 0. Hence we use iff in 
Ea. J18> directly, not the Taylor series. Let us calculate the correction terms in the amplification factor arised from the lens center 
for following cases; a -2,\, and the others. 

Case of or = 2 If the inner density profile is p oc (e.g. the singular isothermal sphere model), the deflection potential is given 
by i^(x) - iffQX (iffQ is a constant) from Ea. J18l l. Inserting this potential t// into Eq.(|5}, we obtain, 

Fiw,y) = I dx'^ exp |^-/ l^yx\ + xj/Q yjxf + x'^ + 6)(l/w)|J , (19) 

where x' = wx. Denoting dFc(w,y) as the leading term of the above integral which is proportional to 1/w, we obtain dF^ as, 

giM.'t('-/2+<A,„Cv)] I 

dF,iw,y) = — forM>|<Aol, (20) 

for \y\ < l^^ol. (21) 



w (/ - (^2)3/2 
„/.rb'2/2+0,„(>O] / 



w (i^l-y^W 

Thus, the contributions to the amplification factor F at the lens center is of the order of 1/w ~ AjM. This is because of the 
singularity in the density profile at the lens center The correction terms dFc in Eas. (l20> and \2\\ represent a diffracted image 
which is formed at the lens center by the diffraction eff'ect. The magnification of this image is of the order of ~ AjM. 



Case of a = 1 If the inner density profile is p oc r ' (e.g. the Navarro Frenk White model), the deflection potential is given by 
i/'(x) - iffQX^ In X {tjjQ is a constant) from Ea.( ll8> . Inserting this ij/ into Eq.(|6}, we have, 

F{w^y) ^ _1 e™[^V2+^„(y)] dx'x'Uyx') e'V(2-K^[i-2^„in(.'/-)]^ (22) 

W Jo 

where x' = wx. By expanding the integrand of Ea. J22> in powers of 1/w, the above equation can be integrated analytically. 
Denoting dFdw, y) as the leading term, we obtain, 

JF,(w,y)=f^e"*'/2^*"<-''". (23) 
Thus, the diffracted image is formed at the lens center with the magnification ~ (A/M)^. 



^ We do not consider the steeper profile a > 2, since the mass at the lens center is infinite. 
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The other cases (0 < or < 1, 1 < or < 2) The deflection potential is ifr{x) - if/o x (i/'o is a constant) from Ea. J18l i. Inserting 
t// into Eq.(|6j, we obtain dFdw, y) as, 

4 (|)"-V„..-3 Hg^, 

where F is the gamma function. 

From the above discussion for the three cases, the diffracted image is always formed at the lens center for the inner density 
profile p oc r^" (0 < a < 2). The magnification of this central image is roughly given by, ^ ~ (A/M)^^". 



5. Results for Specific Lens l\/lodels 

We apply the quasi-geometrical optics approximation to the cases of the simple lens models. We consider the following axially 
symmetric lens models: point mass lens, singular isothermal sphere (SIS) lens, isothermal sphere lens with a finite core, and 
Navarro Frenk White (NFW) lens. We derive the amplification factor F, that in the geometrical optics limit Fgeo, and its correction 
terms dF^ (in §4.1) and dFc (in §4.3) for the above lens models. We define dF for convenience as the sum of dF^ and dFc'. 
dF(w,y) = dF„(w,y) + dFc(w,y). 



5.1. Point Mass Lens 

The surface mass density is expressed as S(^) = M6^(^), where M is the lens mass, and the deflection potential is iff(x) = In x. 
This lens model is most frequently used in wave optics in gravitational lensing of gravitational waves (Nakamura 1998 , Ruffa 
[T999I De PaoUs et al. 12001 2002 Zakharov & Bai-vshev l2002l Takahashi & Nakamura 2003 Yamamoto l2003l Varvella et al. 
120031 Seto l2004t . The amplification factor F of Eq.(|5) is analytically integrated as (Peters il974i Deguchi & Watson [T986> . 



F(w, y) - exp 



nw w 
T ^ '2 



{ln(^)-2M3')}] r(l-^w) iFi(^w,l;^w/), (25) 



where iFi is the confluent hypergeometric function and 0„,(y) - (x^-y)^ /2-lnx„, withx,,, - (y+ yjy'^ + 4)/2. In the geometrical 
optics limit (w » 1) from Eq.(|8} we have, 

Fgeo{^v,y) = - i e"'^\ (26) 

where the magnification of each image is yu± - 1/2 ± (y^ H- 2) / {2y -^y^ -i- 4) and the time delay between double images is 
Ar = y V/ + 4/2 + ln(( V/ + 4 -H y)/( V/ + 4 - y)). In the quasi-geometrical optics approximation, dFc - and dF(- dF„,) is 
given from Ea.( ll4> by, 

3w (x^ + ly{xi - 1) 3w {xi + lyixi - 1) 



where x± - {y + -\/y^ + 4)/2 is the position of each image. The first and second terms in Ea. (l27> are correction terms for the 
magnifications of the two images as discussed in §4. 1 . 

In Fig[na), the amplification factor \F\ is shown as a function of w with a fixed source position y = 0.3. The solid line is 
the result F in Ea.( l25t : the dotted line is the geometrical optics approximation Fg^o in Ea.( l26> : the dashed line is the quasi- 
geometrical optics approximation Fg^o + dF in Eas.( l26> and J27> . For the high frequency limit w » 1, l^^l converge to the 
geometrical optics limit in Ea.(l26>. 

[Fgcof" = Kl + + 2 Ka'-I'/' sin(wAr). (28) 

The first and second terms in above equation ilHi . |//| = |//+| + represent the total magnification in the geometrical optics. The 
third term expresses the interference between the double images. The oscillatory behavior (in Fig^a)) is due to this interference. 
The amplitude and the period of this oscillation are approximately equal to 2|ju+//_|'^^ and wAT in the third term of Eq.j28L 
respectively. For large w (^ 10), these three lines asymptotically converge. 

In Fig^b), the differences between F, Fg^o and Fgeo + dF are shown as a function of w with y - 0.3. The thin solid line is 
\F - Fgeo\, and the thin dashed line is \F - (Fgeo + dF)\. The thick solid (dashed) line is the power of w '(w^^). From this figure, 
for larger w (» 1) F converges to Fgeo with the error of 0(l/w) and converges to Fgeo + dF with the error of 0{l/w^). These 
results are consistent with the analytical calculations in §3 and 4. 
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Fig. 1. (a) The amplification factor \F(w, y)\ for a point mass lens as a function of w with a fixed source position y - 0.3. The solid 
line is the full result F; the dotted line is the geometrical optics approximation dFgeo', the dashed line is the quasi-geometrical 
optics approximation Fgeo + dF. (b) The differences between F, Fgeo and Fgeo + dF for a point mass lens as a function of w with 
y = 0.3. The thin solid line is \F - Fgeo\, and the thin dashed line is \F - {Fgeo + dF)\. The thick solid (dashed) line is the power 
of l/w(l/w2). 



5.2. Singular Isothermal Sphere 

The surface density of the SIS (Singular Isothermal Sphere) is characterized by the velocity dispersion v as, 2(f) - v^/(2f). 
and the deflection potential is il/{x) - x. The SIS model is used for more realistic lens objects than the point mass lens, such as 
galaxies, star clusters and dark halos (Takahashi & Nakamura l2003> . In this model F is numerically computed in Eq.(|6j. In the 
geometrical optics limit w » 1, Fgeo is given by, 

Fgeoiw,y) = - il/Z-l'/^e'-A^ fory < 1, 

= for3;>l, (29) 



where jd± - +1 + 1 /y and Ar = 2y. For y < 1 double images are formed, while for y > I single image is formed. In quasi- 
geometrical optics approximation, dF is given by, 

'^^("'^^ = ^T^2 - i^lT^ + -7rAw2 fory < 1, (30) 

6wy(y+iy 6wy(l-yy w (1 - y^)-^'-^ 

- l^+l'^' +L_L^^n.b^/2^M)} fory>l, (31) 



Swy{y+iy wCy2- 1)3/2 



where <pm(y) = y + 1/2. For y < 1, the first and the second term in Eq.J30> are correction terms for the magnifications of the 
images formed in the geometrical optics (i.e. dF^ in §4.1), and the third term corresponds to the diffracted image at the lens 
center (i.e. dFc in §4.3). For y > 1, the first term of Ea.( l31> is correction term for the magnification dF,,,, and the second term 
corresponds to the diffracted image at the lens center dFc- Thus, in the quasi-geometrical optics approximation, for y < I the 
three images are formed, while for y > I the double images are formed in the SIS model. 

Fig|2]is the same as Fig[2 but for the SIS lens with a source position y = 0.3. In Fig|2ja), the behavior is similar to that in 
the point mass lens (in Fig[na)). The oscillation of \Fgeo\ is because of the interference between the double images, while the 
oscillation of \Fgeo + dF\ is because of that among the three images. As shown in FiglJJb), the errors decrease as w increases, and 
the results are consistent with the analytical calculations. 
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Fig. 2. Same as FigQ] but for SIS lens with a source position y = 0.3. 



5.3. Isothermal Sphere with a Finite Core 

We investigate the effect of a finite core at lens center on the amplification factor We consider the isothermal sphere having a 
finite core. The deflection potential is = (x^ + x^)'^^ where Xc is a dimensionless core radius. In this model, the central core 
of the lens contributes the integral of F in Eq.(|5j. Denoting dFc as the contribution of F at the lens center, we obtain. 



dFc(w,y) 



2mw 



I 



exp 



^ + x;^ + (wxc)^ + 0(i/w) 



(32) 



For wxc $ 1 the above equation is the same as dFc in the SIS model (see Eq.(ll9ll) and dFc oc l/w, but for wxc ^ 1 dFc 
exponentially decreases as w increases. Thus even if the lens has the finite core at the center, the wave does not feel the existence 
of the small core Xc ^ l/w, and the diffracted image is formed similar to the lens without the core. 



5.4. NFW lens 

The NFW profile was proposed from numerical simulations of cold dark matter (CDM) halos by Navarro, Frenk & White J1997I I. 
They showed that the density profile of the dark halos has the "universal" form, p(r) = Ps{f/f.i) \f/rs + l)"^, where r, is a scale 
length and is a characteristic density. The NFW lens is used for lensing by galactic halos and clusters of galaxies. The deflection 
potential is written as (Bartelmann ll996l l. 



2 



arctan 



for X < I, 



for X > 1, 



(33) 



where - 167Tps(DiDis IDsYs is the dimensionless surface density and we adopt the scale radius, not the Einstein radius, as 
the normalization length: ^0 = '"j- The image positions Xj, magnifications jij and time delays Tj are numerically obtained from 
the lens equation, y - x - tfr'{x). In the quasi-geometrical optics approximation, dF is given by. 



dFiw,y)^ - Y Al 



iwT j-mn i 



„/HlvV2+0„,OO] 



(wy2) 



2\2 



(34) 



The first term in Ea. (l34> is the corrections for the magnifications of the images (i.e. dF„ in §4.1). The second term corresponds 
to the diffracted image at the lens center (i.e. dFc in §4.3). 



6. Conclusions 



We studied the gravitational lensing in the quasi-geometrical optics approximation which is the geometrical optics including 
the corrections arising from the effect of the finite wavelength. Theses correction terms can be obtained analytically by the 
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asymptotic expansion of the diffraction integral in powers of the wavelength A. The first term, arising from the short wavelength 
limit /I — > 0, corresponds to the geometrical optics limit. The second term, being of the order of A/M (M is the Schwarzschild 
radius of the lens), is the first correction term arising from the diffraction effect. By analyzing this correction term, we obtain the 
following results: (l)The lensing magnification fi is modified to yu (1 + S), where 6 is of the order of (A/M)^. (2)If the lens has 
cuspy (or singular) density profile at the center p(r) oc r^" (0 < a < 2) the diffracted image is formed at the lens center with the 
mag nification /u ~ (A/M)^-". Thus if we observe this diffracted image by the various wavelength (e.g. the chirp signal), the slope 
a can be determined. 

Acknowledgements. I would like to tiiank Takashi Nakamura, Takesiii Qiiba, and Kazuhiro Yamamoto for useful comments and discussions. 



Appendix A: Numerical Computation for the Amplification Factor 

We present the method for the numerical integration of F discussed in §2. F in Eq.© is rewritten as, 

/r(w,y)= — e'»lv-/2+0»,(y)] dx x e''''^' \ e-™'!-'^^'™^*^'^*^''')!, (A.l) 
27ri Jo Jo 

where 9 is defined as x ■ y = xy cos 9. Changing the integral variable from xio z - x^ 12 in Eq. ( IA.1> . we obtain, 

F(w,y)= fl'z/(z: w,y)e™'^^ (A.2) 
Jo 



/(^ ; vv, y) = — ( d9 ^'"^l . (A.3) 



where the function / is defined by. 

^/.v[j--/2+<J,„(y)l ( 

27r/ 

Thus, we obtain the Fourier integral \A.2l . We present the method for computing Fourier integral in Numerical Recipes (Press et 
g/. n 992, 1. The equation ( IA.2> is rewritten as, 

Jnb poo 
t/z/(z: w,y)e™'+ dzf{z:w,y)e'''\ 
Jb 

dzf(z ■w,y)e : + -—z , (A.4) 

iw [iwY 

where /' - df/db. The first term of the above equation iAAi is evaluated by numerical integration directly. For the series of the 
above equation, we use the integration by parts. The asymptotic expansion in Eq. jA.4> converges for large b. 
Especially for the axially symmetric lens, the function / in Eq. jA.3> is reduced to the simple form as, 

f(z :w,y)^ -we"[>'/2+^".W]yg(^3, ^2^) e'""'. (A.5) 



Appendix B: Asymptotic Expansion of the Amplification Factor for Non-axially Symmetric Lens Models 

We consider the expansion of F in powers of 1 / w for the non-axially symmetric lens model. In this case, the expansion of r(x, y) 
around the j-th image position Xj in Eq.(|9j is rewritten as, 

T(X, y) = + i ^ TabXaXh + ^ ^ dadbdcTiXj, y)XaXtXc + ^ ^ dadbdcddT(Xj, y)XaXbXcXd + O(x^), (B.l) 

a,b a,b,c a,b,c,d 

where Tab - dadtTiXj, y) is a 2 x 2 matrix. We change the variable from Xa to Za - Tjb^abXb in order to diagonalize the matrix 
Tab in the above equation ( IB. 11 1. Here, Aab satisfies, Yia,b TabAacAbd - ^c^cd, where Ac is the eigenvalue of Tab- Using the variable 
z, the equation iB.l\ is rewritten as, 

r(x, y) = Tj + ^ (^Aizl + A2zl) + ^ Mabc ZaZbZc + ^"bcd ZaZbZcZd + 0{z'), (B.2) 

a,b,c a,b,c.d 

where Mabc and Nabcd are defined by, 

Mabc = ^Y^dddcdfT{Xj,y)AadAbeAcf, 

d,e,f 

Nabcd = ^ X ^edfdgdhT(xj,y) AacAbfAcgAdi,. 

e,f,g,h 
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by 
A) 



Doing the same calculations in §4.1 with Ea. (IB.2> . we obtain the amplification factor same as equation il2i . but Aj is given 
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Mf 



+ 3— — — + 3 — , , , ,^ + 



■"^222 



M\^2 



N\\]\ -N\\22 A^2222\ 
+ 2—— + 



1^1 



MA2 \H- 



(B.3) 
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